On the Effects of a Bulk Perturbation on the Ground State of 3D Ising Spin Glasses 
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We compute and analyze couples of ground states of 3D spin glasses before and after applying 
a volume perturbation which adds to the Hamiltonian a repulsion from the true ground state. The 
physical picture based on Replica Symmetry Breaking is in excellent agreement with the observed 
behavior. 
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Studying the scaling behavior of ground states (GS) 
of disordered systems Jlj] turns out to be a powerful tool. 
Here by using a powerful technique based on a bulk modi- 
fication of the Hamiltonian J2|,[| we are able to give strong 
numerical evidence of a Replica Symmetry Breaking like 
behavior of 3D spin glasses. With this study we believe 
we vastly enlarge the scope of where we were chang- 
ing the boundary conditions to analyze the changes in- 
duced in the GS microscopic structure. The open debate 
is about the role of the exact (Replica Symmetry Break- 
ing, RSB) solution of the Mean Field spin glass model, 
where field theory is exact by definition || , in describing 
finite dimensional spin glasses: the validity of a possible 
phenomenological description, the (more or less modi- 
fied) droplet picture should probably be taken as an 
alternative possibility. The RSB picture is characterized 
by the presence of states (for the terminology we refer to 
||) that are locally different one from the other. 

In this respect equilibrium Monte Carlo simulations are 
useful Q, but since it is impossible to thermalize large, 
low T systems, if one is not far enough from T c there 
is the remote possibility j?J that the asymptotic physical 
picture could be disguised because of crossover effects 
(this criticism docs not apply to the off-equilibrium sim- 
ulations |J, and very low T recent simulations also 
claim to disprove it). The remarkable achievement of the 
work of Jj| has been to realize that one can compute GS 
for systems of size comparable to the one of finite T nu- 
merical simulations. The first series of papers has 
been computing couples of GS with periodic boundary 
conditions (PBC) and with anti-periodic boundary con- 
ditions (APBC). A different method, based on the use of 
a bulk perturbation M has been first used by Palassini 
and Young to detect the inconsistency of a (even modi- 
fied) droplet picture: large scale excitations with a finite 
energy cost (first discussed in |l(|) are clearly detected 
by this approach. In the replica approach we expect that 
these large scale excitations correspond to the simulta- 
neous reverse of a large non-compact block of spins: the 
surface to volume ratio of the set of the flipped spins 
should asymptotically go to a non-zero value. In previ- 
ous studies 110,01 this surface to volume ratio was inves- 



tigated: it was found that it was slowly decreasing with 
L and its behavior was compatible with a small negative 
power of L. 

We use here the same method to make stronger the evi- 
dence of the existence of these large scale excitations, and 
we present a detailed quantitative study of their proper- 
ties. We show that they behave in agreement with RSB, 
and that the surface to volume ratio extrapolates to a 
non-zero value in the limit L —> oo. 

We work in three spatial dimensions (D = 3) on simple 
cubic lattices of linear size L, containing N — L 3 spins 
with PBC (we have taken L = 4, 6, 8, 10, 12, 14). We 
consider quenched random couplings J extracted from a 
Gaussian probability distribution with zero expectation 
value and unit variance. We study the behavior of the GS 
after changing the Hamiltonian to include a term which 
repels from the true GS. The Hamiltonian of the first 
system is 



h (W}) 



a(i)J(i, n)o-(i + ji) 



where the sum is over the N sites i and the D directions 
\x (fi is the versor in the /i direction). The Hamiltonian 
of the second system is 

Hi({r}) = H q ({t}) -J^Y, ^o(^o(i+A)r(i)r(i + £) , 



where the o~o(i) are the spins in the GS of Ho. Let 
To(i) be the spins in the GS of H\\ we are interested 
in studying the differences among o~q and To when N 
goes to infinity at fixed e (in this limit the perturba- 
tion is of order 1, while the total Hamiltonian is of order 
N). At this end we define the local site overlap on site 
i as q(i) = (To(i) To(«), and the local link overlap is de- 
fined among nearest neighbor sites (or equivalently on 
links) as qi(i,[i) = q(i) q(i + /i). If two spin configu- 
rations differ by a global spin flip their link overlap is 
identically equal to +1 (while their overlap is equal to 
— 1). We also define the average over the system sites 
and links of the previous quantities: q = A r_1 ^ i g(i), 
qi = {DN)~ 1 J2 ifi qi(i,fi). With these definitions we 
have that Hi({t}) = H ({t}) — eqi. We are consider- 
ing this type of perturbation because the quantity qi is 



1 



sensitive only to local changes of the configuration, i.e. 
it is different from 1 only on the interface among the 
clusters of flipped and non-flipped spins. 

The main effect of the perturbation is to generate ex- 
cited states. It is reasonable to assume that for L — > oo 
the value of the observables computed for fixed value of 
the site overlap are weakly dependent on the value of e. 
This is also suggested by the results of 0], and will be 
discussed in more details in jO]. For this reason we will 
consider the link overlap (and also other quantities) as 
a function of the (square) site overlap, i.e. qi(q 2 ). This 
quantity is defined as the average of the link overlap over 
those systems for which the square site overlaps of the 
two ground state is q 2 . We will indicate the link overlap 
summed over all the systems by q l = J dq 2 qi(q 2 )P(q 2 ), 
where P(q 2 ) is the probability distribution of q 2 . We 
have computed 10 4 GS at L = 4 and 6, 2000 at L = 8, 
1500 at L = 10, 500 at L = 12 and 14, with e = 4^6, 
using a genetic algorithm which will be described in [fll| . 

Given the importance of the issue, we start by dis- 
cussing the dependence of qi{q 2 ) over q and L. If we 
plot qi(q 2 ) versus q 2 , we notice that qi{q 2 ) and q 2 are 
very strongly correlated. In order to clarify this effect 
we define the quantity (for all systems such that q ^ 1) 
R L {q 2 ) = (1 - qi(q 2 ))/(l - q 2 ). For fixed L R L (q 2 ) is 
nearly q 2 independent in a large region where q is not 
too large (only close to q = 1 the dependence becomes 
slightly more noticeable). The q 2 dependence of R(q) in 
this large q region becomes milder with increasing L. We 
fit R,L(q 2 ) to the form A(L) + B(L) q 2 using the data in 
the range q 2 £ (0,0.3): the fits are very good, and the 
extrapolations are all very smooth. B(L) is in all cases 
a few percent of A(L). In fig. [I] the upper points are for 
A(L) = R L (0) = 1 - qi{0) versus L~ l . Both a linear fit 
(which extrapolates to a non zero value) and a power fit 
to the form CL~ a , where a turns out to be 0.44, follow 
the data in a qualitative way, but they are unacceptable 
because they lead to a very high value of x 2 (i- c - to a 
normalized % 2 = O(20)). In fig. [I] we plot a second or- 
der polynomial fit that is very good (it has a normalized 
X 2 ^ 0.6) and leads to a value of A(L = oo) = 0.30±0.01. 
We also plot a fit to the form A(L) = a + bL~ x , that is 
also good (it has a normalized x 2 — 0.3) and leads to a 
value of of A(L = oo) = 0.24 ± 0.01 and A = 0.72 ± 0.02. 
The difference with the 0.30 obtained by the polynomial 
fit gives an idea of the size of the systematic error over 
the extrapolation. 

The values extrapolated to L = oo are an estimate of 
1 — qi(0), and they have to be compared with the value 
0.245 ± 0.015 that we have estimated in Q: the agree- 
ment is good, even if here the systematic error is large. 
We can conclude that there is evidence that 1 — qi (0) does 
not extrapolate to zero in the infinite volume limit and 
the dependence of qi(q 2 ) is well approximated (for q 2 not 
too close to 1) by the form qi(q 2 ) = 1 — A(l — q 2 ). The 
previous formula holds in the infinite range SK model, 



with A = 1. The fact that in 31? we find a low value 
(A ~ 0.3) suggests that A goes to zero when one ap- 
proaches the lower critical dimension. 

The lower points and curves of fig. [j] are for the ex- 
pectation value of the link overlap averaged over all 
q 2 values. The curves are again a polynomial and a 
power law best fits, that here extrapolate to a lower value 
(0.07 ± 0.04 for the power law and 0.155 ± 0.015 for the 
polynomial fit). The error is larger than in A(L) (of a 
factor close to 2) but apart from a shift (since 7ji ^ qi(0)) 
the two curves are similar. 
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FIG. 1. The value of A(L) (upper points) and of q l as a 
function of L _1 . The curves are for a power fit and a polyno- 
mial fit. 

Let us stress that A(L) is important also since qi(q 2 ) 
is not expected to depend on e when L — > oo. The e- 
perturbation allows to explore low energy states with an 
energy gap of order 1. The function P(q 2 ) changes with 
e p2] and consequently also the value 7ji changes: still 
the functional relation that relates qi to q 2 is unchanged. 

It is also interesting to look at the variance squared 
of R L {q 2 ), i.e. to a 2 R)L (q 2 ) = (R 2 L (q 2 )) - (R L (q 2 )) 2 . 
We proceed as for the the best fit of RL{q 2 )- First we 
fit <y 2 R ^ L {q 2 ) to the form Aw{L) + R>w{L) q 2 , using the 
data with q 2 < 0.3: the extrapolations are very good. 
Also here for a given value of L we find that cf 2 R ^ L {q 2 ) 
is nearly q 2 independent apart from the region close to 
q 2 = 1. Also the q 2 dependence of ^^^(q 2 ) in this re- 
gion becomes weaker mild with increasing L. 

We plot Aw{L) versus 1/L in fig. ^. Together with 
the data we show a best fit to the form A\y(L) ~ C L~ 5 , 
that has a low x 2 > and 8 — 1.35±0.03. The fact that a fit 
to a zero asymptotic value is very good shows that qi as 
a function of q does not fluctuate as L — > oo : this is the 
same situation that we find in the mean field RSB theory, 
with a more complex dependence of qi(q 2 ) over q 2 . The 
asymptotically vanishing value of Aw(L) explains why 
the errors on qi(q 2 ) are smaller that those on ql. 

We discuss next the overlap-overlap correlation func- 
tions. The difference among a droplet and a RSB scenario 
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becomes very clear if we consider the overlap-overlap cor- 
relation function at fixed value of q 2 in a system of side L, 
i.e. the correlation function obtained by selecting those 
systems which have a ground state with a fixed value of 
the mutual overlap: 
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FIG. 2. The value of A W (L) as function of L' 1 . The curve 
is a power fit to a zero asymptotic value. 

C{x,L\q 2 ) = {q{x)q{Q)) {<f) . 

In the RSB case we expect that for large L and fixed x 
C(x,L\q 2 ) — > G{x\q 2 ) + q 2 , where the function G(a;|g 2 ) 
depends on q 2 and for large x should behave as G(x\q 2 ) oc 
x^^ q K (f){q 2 ) may depend on q 2 , although there are in- 
dications that it takes at most only two different values, 
one at q = and one at q ^ jl3| |. In the region where 
x and L are large we expect that 

C(x, L\q 2 ) » g{r\q 2 ){- + -J—)^ 2 ) + q 2 , (1) 
X 1j — X 

where r = xL~ l . We have used the factor (- + jtt" ) 
instead of the simpler i in order to implement the sym- 
metry implied by the presence of PBC. 

The situation described by a droplet picture is differ- 
ent: in this case, for large L and fixed x, C(x, L\q 2 ) — > 1. 
One can argue that in the droplet model one should 
have 0(q 2 ) = and g(0\q 2 ) + q 2 ~ 1, i.e C{x,L\q 2 ) » 
D(xL~ l \q 2 ) with D(0\q 2 ) = 1. 

We will discuss our results for small values of q: we 
consider the g range q < q ma2; =0.1. The results are not 
very sensitive to the exact position of the q ma x cut (of 
course if q m ax is too small the statistics fades) . In fig. ^ 
we present our data for g(r) as a function of w = r(l — r) 
for = 0.4 (we prefer to plot g versus w since because 
of the PBC both g and w are left invariant when r goes 
in 1 — r) . The value of q 2 we use to subtract the discon- 
nected part is the average of q 2 in the interval q < q max . 
We have selected <fi such to optimize the quality of the 
scaling behavior. The line is for the best fit to a second 
order polynomial. 

The scaling is very good, and the estimate of <f> — 0.4 
reliable. A droplet-like scaling is excluded by our data. 



We notice that C(x, L\0)\ x= i = qi(0) and the scaling form 
(|lf) accounts for the L-dependence of qi(0). If we assume 
that (|]) is exact also for x = 1 (a too strong assumption) 
we get that C(l,oo|0) = <?(0|0). Using the quadratic fit 
in fig. U to extrapolate the data at r — we get in the 
infinite volume limit 1 — qi(0) = .33 ± .02, which is in 
reasonable agreement with the previous estimates, given 
the crudeness of the approximation. 
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FIG. 3. The overlap-overlap correlation function minus q 2 
divided times the factor f i + jj^J (i.e. the function g de- 
fined in (Q)) as a function oi w = j- — (r) 2 - The line is for 
the best fit to a second order polynomial. 

When we analyze the correlation function of the link 
overlap we also select events at a fixed value of q 2 : as 
in the previous case we will consider a cut of the form 
M < qmax = -1. We only measure a simple correlation 
function of 1 — qi(i) (we only look at correlations in a 
given plane, and average the two contributions coming 
from separating in the direction 2 links in the direction 
1 and in the direction 1 links in the direction 2 Q), that 
we call Q(x,L), which we plot in fig. |] for various L 
(x > 1). The different sets of points are plotted ver- 
sus s(x,L) = (1/x + 1/(L — x)) x . If we chose the value 
A = 2.1 the data at given L can be very well fitted by a 
linear functions, i.e. the lines in fig. ^: 

Q(x, L) — a{L) + b(L) s(x, L) . (2) 

We plot in fig. || a(L) and b(L) from (||) and the best 
linear fits. If for L — > oo 1 — <?z(0) would go to zero, this 
correlation function should become identically zero and 
consequently the two functions, a(L) and b(L), should 
extrapolate to zero when L — > oo. However the slope 
b(L) does not show any tendency to go to zero: the 
fact that it increases (slightly) with L is a substantial 
indication that it does not converge to zero in the in- 
finite volume limit. The intercept is fitted very well 
by a 1/L behavior without corrections and it extrapo- 
lates to a non zero value. Moreover the correlation func- 
tion at large distances should be the square of the ex- 
pectation values and consequently we have the relation 
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ci(L)\l =00 = (1 — qi(0)) 2 . This relation gives a value of 
1 — qi (0) = 0.26 ± 0.02, that compares very nicely to the 
values we have discussed earlier. 

Let us start summarizing our results. We have applied 
a bulk perturbation, and shown that the results substan- 
tiate in a perfect way the ones obtained by changing the 
boundary conditions 0]. The quantity 1 — qi(0) has a 
non-zero value that compares very well to the one deter- 
mined in 0. We have also determined that qi behaves 
very similarly, and we have estimated its value in the 
infinite volume limit (four different estimates, affected 
by different systematic errors are .24 ± .01, .26 ± .02, 
.30 ± .01, .33 ± .02, the errors being purely statistical). 
We have shown that very plausibly in the limit L — *■ oo 
qi is a well behaved function of q 2 , in strict analogy to 
what happens in the mean field RSB theory. We have 
analyzed site overlap and link overlap spatial correlation 
functions, and we have shown they perfectly obey scal- 
ing laws implied by a RSB scenario. Again, this results 
is quantitative, in which it leads to an estimate of qi(0) 
coherent with the ones given before. 
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FIG. 5. a(L) (upper points) and b(L) (lower points) from 
equation (0) and the two best linear fits. 





0.6 r 




0.55 




0.5 




0.45 


A 






0.4 


5" 






0.35 


"Si 




V 


0.3 




0.25 




0.2 




0.15 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

(l/d+l/(L-d))**2.1 

FIG. 4. Points are for the link overlap auto correlation 
functions versus w 2 ' 1 . Lattice sizes are from L = 4 (upper 
points and line) to L = 14 (lower points and line). Lines are 
for linear best fits. 

The RSB scenario is perfectly compatible with the 
whole set of data. These results are consistent with the 
recent detection of large scale excitations in spin glass 
ground states jlC|fl. The behavior we find connects 
smoothly with the results obtained by simulations at fi- 
nite temperature || (e.g. the overlap-overlap correlation 
decays as the distance to a power (f> — .4, which is not far 
from the value obtained from simulation at finite temper- 
ature). We think that we have conclusively shown that 
the ground state structure of 3D Ising spin glasses with 
Gaussian quenched random couplings is not trivial. 
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